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ABSTRACT 

• The current methods available to estimate gravitational shear from astronomical images of 

galaxies introduce systematic errors which can affect the accuracy of weak lensing cosmolog- 
Q , ical constraints. We study the impact of KSB shape measurement bias on the cosmological 

i—{ • interpretation of tomographic two-point weak lensing shear statistics. 

C/3 I We use a set of realistic image simulations produced by the STEP collaboration to derive 

^ ■ shape measurement bias as a function of redshift. We define biased two-point weak lensing 

statistics and perform a likelihood analysis for two fiducial surveys. We present a derivation 
of the covariance matrix for tomography in real space and a fitting formula to calibrate it for 
non-Gaussianity. 

We find the biased aperture mass dispersion is reduced by ^ 20% at redshift ^ 1, and 
OO . has a shallower scaling with redshift. This effect, if ignored in data analyses, biases erg and wq 

OO ' estimates by a few percent. The power of tomography is significantly reduced when marginal- 

ising over a range of realistic shape measurement biases. For a CFHTLS-Wide-like survey, 
. [r^m, cs] confidence regions are degraded by a factor of 2, whereas for a KIDS-like survey the 

factor is 3.5. Our results are strictly valid only for KSB methods but they demonstrate the need 
OO ', to marginalise over a redshift-dependent shape measurement bias in all future cosmological 

■ analyses. 
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1 INTRODUCTION 

Weak gravitational lensing by large-scale structure, or cosmic shear 
(see Munshi et al. 2008 for a review), is a powerful tool to inves- 
tigate dark energy and the large-scale distribution of dark matter 
dAlbrecht et al. 20061 [Peacock et aO OOe i. In the last decade, two- 
point shear statistics have been successfully measured and used to 
constrain parameters of the matter power spectrum (from the early 
constraints of van Waerbeke et al. 2001 to the recent results of Fu et 
al. 2008) and dark energy jjarvis etal."2006 . Hoekstra et al. 20061 
ISemboloni et al. 2006IISchimd et al. 2007l|Kilbinger et al. 2008) . 

The information contained in the cosmic shear signal is more 
efficiently exploited if the redshift of the source galaxies is avail- 
able. In particular, the measurement of shear statistics in red- 
shift bins, or tomographic weak lensing, can greatly improve 
cosmological constraints dHu 19991 , especially for dark energy 
dTakada & Jain 2004I I. since the weak lensing signal scales differ- 
ently with redshift for different models of dark energy evolution. 
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The accuracy of cosmic shear constraints is affected by sys- 
tematic effects. Systematics arise primarily from intrinsic align- 
ments, redshifts uncertainties, and shear measurement errors. In- 
trinsic alignments are non-cosmological sources of shear correla- 
tion, caused by an intrinsic correlation between the orientation of 
galaxies dCrittenden et al. 200 11 1 or by shear-shape correlation, i.e., 
by a correlation between galaxies shapes and surrounding density 
fields ( [Hirata & Seljak 2004| l. These effects have been modeled, 
eg., in Heymans et al. (2006b) and have also been constrained with 
data (Mandelbaum et al. 2006V In particular, it has been shown 
(Hirata et al. 2007 ; Bridle & Ki ng 2007[ l that ignoring the effect of 
the shear-shape correlation on the amplitude of the two-point shear 
statistics heavily biases the constraints on the equation of state of 
dark energy. The impact of redshifts uncertainties on cosmic shear 
has also been analysed in several studies (Ish ak & Hirata 20051 
|Ma, Hu & Huterer 2006 , van Waerbeke et al. 20(56) 

This paper deals with the third primary effect: shear mea- 
surement errors due to the lack of an unbiased method to estimate 
gravitational shear from astronomical images of galaxies. It aims 
to derive the effect of shear measurement bias on dark energy 
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constraints both with and without tomography. The estimation 
of gravitational shear is technically challenging and the methods 
available today do not seem to do better than the percent level of ac- 
curacy. The Shear Testing Programme (STEP), and more recently 
the gravitational lensing estimation accuracy test (GREAT08) 
(Bridle et al. 2008), represent to date the largest collaborations 
aimed at testing and improving the accuracy of the existing shear 
measurement methods with the use of realistic sets of simu- 
lated images. Current methods are based on the KSB approach 
( [Kaiser, Squires & Broadhurst 1995[ [Luppino & Kaiser 1997[l, 
shapelet decomposition dBernstein & Jarvis 20021 

[Refregier & Bacon 2003[ | Kuijken 2006[ >, or Bayesian tech- 
niques I Kitching et al. 2008b i. The STEP collaboration has shown 
that the methods developed in the past years underestimate on 
average the shear signal by a few percent ( [Heymans et al. 2006a [ 
[Massey et al. 2007] l, implying the measurement of two-point shear 
statistics is underestimated by two to ten percent. Moreover, 
the STEP collaboration showed that this bias is not constant but 
depends on the brightness and size of the galaxies. Indeed, for most 
methods the measured shear is overestimated for bright galaxies 
and greatly underestimated for faint galaxies. This may be a 
limitation of the accuracy with which one can derive cosmological 
constraints using cosmic shear statistics in future ground-based 
and space-based astronomical surveys such as KIDS, LSST, 
Pan-STARRS, DES, JDEM or EUCLID. 

There are in the literature various studies dealing with the 
impact of shear measurement errors and PSF modeling on cosmic 
shear measurements i jHirata & Seljak 2003[ [Hoekstra 2004. 
Ivan W aerbeke et al. 2005 ; [Jain, Jarvis & Bernstein 2006] 

IPaulin-Henriksson et al. 2008b . There are also several cosmic 
shear forecasts of cosmological parameters which include and 
model generic types of systematic errors, such as multiplicative 
and additive errors in shear measurements tHuterer et al. 2006. 
[Amara & Refregier 2007| [Kitching, Taylor & Heavens 2008a^ . In 
these previous works the bias has been assumed to be a generic 
function of the redshift; in this paper we adopt a different approach, 
consisting in explicitly quantifying the bias as a function of redshift 
using realistic simulations, as opposed to using generic modeling. 
In order to derive the expression of the bias as a function of redshift 
we reanalyse the STEP2 simulations using a KSB pipeline. We 
observe that the bias depends on the characteristics of the galaxies 
such as magnitude and size and we use these dependencies to 
derive a realistic expression of the bias as a function of redshift. 
In order to evaluate the impact of this bias on the estimation of 
cosmological parameters we perform a likelihood analysis of to- 
mographic two-point cosmic shear in a space of cosmological and 
bias parameters, including small-scale non-Gaussian corrections 
in the covariance matrix. We note that even though we model the 
bias using a particular implementation of KSB, the results we find 
are general as all KSB methods tested by the STEP collaboration 
show the same behavior. 

The paper is organized as follows: in Section |2| we describe 
the results of the STEP2 analysis and the procedure to estimate the 
bias as a function of redshift. Section [3] defines the tomographic 
two-point shear statistics in the presence of a redshift-dependent 
bias. In Section |4| we define the space of cosmological and shear 
measurement bias parameters where the likelihood analysis will be 
performed. Section|5]presents the correction for non-Gaussianity to 
be applied to the covariance matrix of tomographic two-point cor- 
relation functions computed for a Gaussian shear field. The deriva- 
tion of this covariance in the real space is detailed in the Appendix. 



The results of the likelihood analysis are presented and discussed 
in Section|6|and we conclude in Section|7| 



2 STEP2 ANALYSIS 

In the last decade a large effort has been made in order to estab- 
lish the accuracy with which one can estimate shear from images 
of galaxies. The Shear Testing Programme is a collaboration which 
aims to test and improve methods for PSF correction using realis- 
tic sets of simulations. The STEPl and STEP2 simulations have a 
depth similar to the images of the CFHTLS-Wide survey, which is 
to date the largest deep weak lensing survey available. The results 
of the analyses on the simulations can thus be immediately used to 
assess the accuracy on the cosmological constraints from real sur- 
veys. The galaxies in the STEP2 set of simulations have a bright- 
ness profile built using a shapelets decomposition of galaxies ob- 
served with the Hubble Space Telescope j Massey et al. 2004) ; the 
resulting profiles are more realistic than the STEPl simulations in 
which elliptical galaxies are characterised by de Vaucouleurs pro- 
files and spiral galaxies are characterised by a bulge and a projected 
disk. 

We recall here the main characteristics of the STEP2 simula- 
tions referring the reader to Massey et al. (2007) for further infor- 
mation. The STEP2 set of simulations is composed of six subsets, 
each characterised by a different PSF. Each subset is identified by 
an alphabetical letter from A to F. The PSFs A B and C are typ- 
ical Subaru PSFs. The PSF D and E are highly elliptical, aligned 
along the x and y axis, respectively. The PSF F is a circularly sym- 
metric PSF. Furthermore the PSF is constant across the field; the 
seeing size is ~ 0.6 arcsec, apart from PSF C for which the seeing 
is ~ 0.8 arcsec. For each PSF there are 64 pairs of 7' x 7' images 
with pixel size of 0.206 arcsec. Each pair corresponds to one image 
and its rotation by 90 degrees. A shear field and a PSF are applied 
to the image pairs. The image pairs have been generated so that one 
can estimate the shear using a galaxy and its rotated pair, reducing 
the effect of shape noise on the measured shear. The STEP2 sim- 
ulations are therefore more suitable for our study than the STEPl 
simulations both because the galaxy profiles are more realistic and 
the bias can be determined to a higher level of accuracy. For this 
reason we will consider the STEP2 set of simulations as our main 
dataset, but we will briefly discuss the results we obtain using the 
STEPl set of simulations at the end of this section. 

We analyse the STEP2 simulations using the KSB implemen- 
tation used in Fu et al. 2008 but with the following significant 
changes: 

• we do not introduce a weighting scheme, i.e. all the objects 
have the same weight 

• we use the trace of the shear polarizability tensor Pg in the 
shear estimator for each object 

• we select objects only according to their signal-to-noise ratio. 
The signal-to-noise is defined by the v parameter provided by the 
imcatQ shape measurement software. 

We introduce these modifications as we want to keep the PSF cor- 
rection method as basic as possible avoiding arbitrary "ad hoc" 
choices which might change the results significantly. For each PSF, 
we analyse the 64 pairs of images, merging the catalogues of the 
image pairs so that the final catalogue is free of spurious detections. 



| http : / / www ■ if a ■ hawaii ■ edu/~kaiser developed by Nick 
Kaiser. 
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Figure 1. Left panel: values of the multiplicative bias mi as a function of magnitude averaged on the 64 image pairs. Error bars represent the standard 
deviations around the average value. Different colors refer to different PSF sets. The dashed lines show the values of the multiplicative bias factor mi 
averaged on the whole catalogue. Right panel: same as the left panel but for the second component m2 . 
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Figure 2. Left panel: average value of the multiplicative bias factor m as a function of the redshift for each set of PSF. The value has been obtained by 
averaging over the one hundred bootstrap realisations built as indicated in the text. Error bars represent standard deviations of the one hundred realisations. 
The dashed lines represent the best hnear fit to the points. Right panel: same as the left panel for the additive bias constant c. 



As suggested in Massey et al. (2007) we describe the differ- 
ence between the measured shear ■y and the input shear signal "f^,.^^ 
using the multiphcative m — (mi, 7712) and additive c — (ci, C2) 
bias factors: 

7i (^) - 7«,truc(0) = rrii 7i,truc(0) + Cj i = l,2. (1) 

The estimated shear at a position 9 is affected by a multiplicative 
bias m. The main reason for this calibration error is that KSB meth- 
ods do not completely correct for the isotropic kernel of the PSF, 
leaving thus a residual which is seeing and magnitude dependent. 
Since the seeing in the STEP2 simulations is constant across the 
image, m is independent of position. The additive bias c comes 
from incomplete correction of the PSF anisotropy. This originates 
from the approximative nature of the KSB method and also from 
the residuals due to the modeling of the variation of the PSF across 
the image, which depend on the distribution of the observed stars. 



The obtained m and c biases depend on the magnitude and on the 
size of the galaxies. We define the size of an object as the charac- 
teristic size r;i of the Gaussian profile which best fit the brightness 
profile as defined by the IMCAT software. We verify that our results 
are similar to those obtained by STEP2 by measuring the values of 
m and c as a function of magnitude for the two ellipticity com- 
ponents (see Fig.[T). For example, the value of the m components 
averaged over all galaxies lies between —7% and +2% depending 
on the PSF subset, with a mean over PSFs of ^ —3%. 

We want now to convert the bias as a function of magnitude 
and size into a bias as a function of redshift. To do so, we need 
to find the average redshift distribution of the galaxies belonging to 
each size and magnitude bin, using the CFHTLS-Deep photometric 
redshift catalogue dllbert et al. 20061 1. For this we first convert the 
STEP2 magnitudes into the corresponding CFHTLS-Deep values. 
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Figure 3. Left panel: size distribution of the galaxies in each of the STEP2 sets (solid histograms). The distribution of sizes of the CFHTLS-Wide galaxies 
is shown for comparison (dashed histogram). The size of the galaxies in the CFHTLS Deep catalogue have been converted into the STEP2 simulations pixel 
units of 0.206 arcsec. Right panel: distribution of magnitudes for each of the STEP2 PSF sets (solid) and for the CFHTLS-Wide galaxies (dashed). 



Table 1. Best- fit values of the redshift bias parameters am , bm ,o,c,bc and 
their relative errors for each PSF set of the STEP2 simulations. These values 
are used to plot the lines shown in Figure|2] The last Une shows the averages 
over all PSFs. 



PSF 


am(10-2) 


bm (10-1) 


ac (10-4) 


be (10-4) 


A 


-2.5 ±0.7 


-0.10 ±0.06 


-3.5 ± 1.4 


-9.8 ± 1.3 


B 


-1.8 ± 1.1 


0.28 ±0.09 


-0.4 ± 1.0 


-13.1 ± 1.2 


C 


-1.0 ± 1.0 


-0.23 ±0.08 


-5.2 ±2.0 


-6.0 ± 1.7 


D 


-1.4 ± 1.0 


-0.41 ±0.07 


-11.1 ± 2.7 


-52.6 ±2.2 


E 


-1.7 ±0.8 


-0.41 ±0.06 


-2.5 ±2.3 


16.7 ± 1.9 


F 


-0.9 ±0.7 


-0.11 ±0.06 


-2.3 ± 1.7 


-0.8 ± 1.3 


mean 


-1.6 ±0.4 


-0.02 ±0.03 


-4.2 ±0.8 


-10.9 ±0.7 



This is done by simply adding the photometric zero-point specified 
in the header of the images. We then find the conversion between 
the observed size of an object in each set of simulations and the 
size it would have in the CFHTLS-Deep catalogues. A conversion 
is needed because the pre-seeing size is the reference size which 
one should use to compare populations of galaxies between cata- 
logues with different seeing. This is in principle not straightforward 
because the observed size depends on the seeing through a convo- 
lution. We check the difference of the seeing values comparing the 
locus of the stars in the simulations and in the CFHTLS-Deep cat- 
alogues. To do this we express the size of the stars in the CFHTLS 
Deep catalogue in units of pixels of the STEP2 simulations (the 
pixel size of the STEP2 simulations is 0.206 arcsec whereas the 
pixel size of the CFHTLS images is 0.186 arcsec). We find that 
the difference between the star size in each set of simulations and 
the star size in the CFHTLS-Deep catalogue is always smaller than 
1/5 pixels except for the PSF C where the difference is about half 
a pixel. From this, we conclude that we do not need to reassign the 
size of the STEP2 galaxies. 

We chose a magnitude range 21.5 < mag < 24.5, for which 
a survey with the same depth of the simulations is complete. We 
divide the galaxies in bins of size and magnitude. By the defini- 
tion (Eq.[T), m and c are two-dimensional vectors but for our pur- 
poses we set both components equal to the average between the 
two components and define scalar quantities: m = ^{mi + m.2) 
and c = ^(ci + C2). For each bin, i, containing rii galaxies, we 



compute the average bias parameters rrii and Ci and respective un- 
certainties Srrii and Sci. We divide the CFHTLS-Deep catalogue 
into the same magnitude-size bins. Each galaxy in the simulation is 
assigned a value nij and Cj , extracted from two Gaussians of mean 
m,; and Ci and dispersion SrUi and Sci, and a redshift Zj extracted 
randomly from the CFHTLS-Deep galaxies belonging to the same 
bin. We bin the result in redshift bins Zk, and the total bias in each 
redshift bin zt is given by: 



En [. 



c{zk) 



(2) 



where is the number of galaxies in the redshift bin. We repeat 
this procedure one hundred times. The final value of m(z) and c{z) 
and the relative errors are obtained averaging over the one hundred 
realisations. 

We chose to fit m{z) and c{z) using a linear function: 



m{z) 



z + b„ 



c{z) 



■ z + bc 



(3) 



We perform the fit using only objects with redshifts z < 2.5. We 
show the values of the multiplicative and additive biases as a func- 
tion of the redshift in Figure |2] where the best linear fit to the data 
is also shown. The best-fit parameter values and relative errors are 
summarized in Table [T] Overall the shear is generally underesti- 
mated by a few percent. The multiplicative bias parameter 771(2) 
has a negative slope am which is roughly 1-2% for all the PSF 
types. The additive constant c{z) shows a slope which is consistent 
with zero except for PSF D. The value of the constant be is also 
significant in particular for PSF D and E. We repeated the analysis 
assuming the bias to depend only on the magnitude, i.e., dividing 
the catalogue only in bins of magnitude, and found similar results 
for the functions 777(2) and c(z). This is due to the fact that the size 
distribution as a function of redshift does not change significantly, 
so even though the bias slightly depends on the size the average 
bias is roughly the same for each redshift bin. 

One should expect the values of the bias parameters to change 
if the selection criteria change. For example, including low signal- 
to-noise objects increases the overall bias and also affects the de- 
pendence on redshift, a„i. In order to explore pessimistic cases we 
removed the signal-to-noise threshold (spurious detections are still 
removed by merging each object with its rotated companion). The 
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biases increase significantly. The mean m over all PSFs is now 
larger than 10%, whereas the redshift-dependence has a slope a„i 
between —6% and —4%, depending on the PSF. This is a conse- 
quence of the fact that lowering the signal-to-noise threshold we in- 
clude in the final catalogues a larger number of faint objects which 
are preferentially at high redshifts. This result shows that when us- 
ing simulations to find a recalibrating factor, capable of compensat- 
ing an average shear measurement bias, the factor obtained can only 
be safely applied to a similar population of galaxies. Likewise, our 
results are only strictly valid for surveys with size and magnitude 
distributions similar to the STEP2 catalogue. As we already men- 
tioned, this is likely to be the case for the CFHTLS-Wide, which 
has a similar depth to the STEP2 simulations. We show in Figure 
[3]that the size and magnitude distributions of STEP2 indeed match 
those of the CFHTLS-Wide. 

We analysed the STEPl sets of simulations and derived their 
redshift-dependent biases. We notice that the distribution of galax- 
ies in STEPl includes a higher fraction of large galaxies than real 
data for the same depth such as the CFHTLS-Wide data. We find 
that the average multiplicative bias is ^ 2%, averaged over the 
PSFs subsets with much less dispersion between PSFs than for 
STEP2. It is characterised by a slope value a™, always smaller 
than 1%, indicating that the accuracy of the shape measurement 
degrades less rapidly when the signal-to-noise ratio decreases in 
comparison to the STEP2 simulations. 

The analysis of STEPl and STEP2 simulations give the same 
qualitative results: the average bias is negative and the bias in- 
creases with redshift. The amplitude of the bias measured in the 
STEPl simulations is smaller than the bias in the STEP2 simula- 
tions. The difference in the results is not surprising as the shape 
and the amplitude of the measured bias depends on the character- 
istics of the simulations. We conclude that in order to quantify the 
effect of the bias on real data one should use simulations that are 
as realistic as possible and as we already pointed out the STEP2 
simulations are more realistic than the STEPl simulations. 



3 GRAVITATIONAL LENSING STATISTICS 

In the next sections we explore the impact of the shear bias on cos- 
mological constraints obtained with two-point shear statistics, both 
with and without considering redshift binning. To do so we de- 
rive an expression for two-point shear statistics which include the 
redshift-dependent bias of the shape measurement. 

From Eq. l[T), the measured correlation function is related to 
the true correlation function via 

(7 (9) 7*(e')) = ([1 + m{z)][l + m(z')]7truc (6) ■y:,uc{0')) + 
+ {c{z)c{z')}. (4) 

Since the two types of biases, m and c, are independent there are 
no mixed terms in Eq. The first term of the right-hand side 
is identical to the correlation function, (Kmod(6') Kmod(6''))' °f ^ 
modified convergence Kmod = (1 + m)Ktrue. The power spectrum 
of the modified convergence relates to the matter power spectrum. 
Using the Limber approximation jBartelmann & Schneider 2001b . 
the relation is: 



\fK{w) ) 



with 

gi[w) = / dw [l + m{w (z))]p^,i{w ) — , , — . 

(6) 

is the comoving angular-diameter distance and m{w) is the 
multiplicative bias written as a function of the comoving distance 
w. Eqs. ^ and (|6j take into account tomography, describing auto- 
and cross-correlations between redshift bins [zinf.i, ^sup.i] where 
the distribution of the sources Pw,i{'w) is zero for w < «;(zinf,i) 
and w > ^(^sup.i). Notice the effect of the multiplicative redshift- 
dependent bias is equivalent to a change in the source redshift dis- 
tribution. Since 1 + m{z) is a decreasing function (cf. Fig.|2]l, the 
effective redshift distribution is shallower. 

The modified two-point shear measured in apertures can be 
computed as usual, by integrating the shear-shear correlation func- 
tion, Eq. l[T), with the appropriate filters. For the top-hat variance at 
an angular scale we obtain 
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(7) 



The first term integrates the modified convergence power spec- 
trum, Eq. using the top-hat filter in Fourier space, W\ (r\) = 
Ji {ri)/{nri), where Ji {rj) is the Bessel function of first order. The 
second term of Eq. l|7j is the contribution from the additive bias. It 
is the integral in the aperture of the second term of Eq. (Til, using 
the top-hat filter in the real space. Notice that due to the redshift 
dependence of c the integral over the aperture also includes an inte- 
gration in the radial direction, weighted by the source distribution. 
If the correlation (c (z) c {z') ) is constant across the image, the ef- 
fect of the additive redshift-dependent bias is to add a constant to 
the two-point shear statistics. 

Similarly, the variance of the aperture mass is 



{M^^m = 2n / dssP4s)[W2{s&)Y 
Jo 



(8) 



with W2{'i]) — 12 Jiiji) / {-Krf), where J4,{r}) is the Bessel func- 
tion of fourth order. In this case there is no contribution from the 
additive bias because the filter of the aperture mass dispersion in 
real space is a compensated filter and thus insensitive to a constant 
shear. 

Figure |4] shows biased and unbiased two-point shear statistics 
for three redshift bins. Only the three auto-correlations are shown. 
For this example we plot Eqs. ([TJ and ^ using a cosmological 
fiducial model, and bias parameters similar to PSF D and E, which 
are the most extreme cases. We parameterise the distribution of the 
sources, pn,{w), as in Benjamin et al. (2007) 



p{z) 



N z!^ + 20 



(9) 



(5) 



normahzed to I, with a = 0.73, fi = 4.52, 20 = 0.80, cor- 
responding to a median redshift z„i = 0.77. These values were 
derived by fitting the density of galaxies nu found in Section|2]for 
each PSF, using only galaxies with z < 2.5. We find similar results 
for all PSFs. 

The left panel of Fig. |4] shows that the modified aperture mass 
variance decreases by a scale independent factor of 15% to 30%, 
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Figure 4. Biased and unbiased two-point shear statistics as function of angular scale i9 for three redshift bins. The unbiased signal (solid lines) is compared 
with the signal that would be measured on a survey with very elliptical PSF (dashed lines) characterised by the following bias parameters: am = —0.03, 
bm = -0.06, ac = -0.001, be = -0.006. Left panel: for {Mfp(i9)). Right panel: for the top-hat variance (7^(1?)) . The extra lines (dash-dotted) show a 
case where the additive bias only correlates on small scales. See text for details. 



increasing with redshift. This increase produces a shallower scaling 
of the amplitude of {M^p{9, z)) with the redshift. Modeling the 



redshift scaling with a power law we find it changes as: 

{Mipi^)) (X z^-°^ (no bias) 

(Mfp(i?)) oc {am = -0.03, 6„ = -0.06). 

The right panel of Fig. |4] shows that the top-hat variance is 
dominated by the additive constant. Indeed, assuming the additive 
bias is constant across the image and correlates at all scales, the 
second term of the right-hand side of Eq. Q writes 



C 



/ 



li>(ZBup) 



dw c {'w{z)) Pw{w) 



(10) 



"(Zinf) 

A constant bias of |6cl ~ 10"'', which is the average value be- 
tween the PSFs shown in the Table [T] yields C ~ 10~^ which is 
significant when compared to the cosmic shear signal, especially 
at the lower redshift bins, invalidating any cosmological interpreta- 
tion. This is in agreement with Huterer et al. (2006) where it was 
found that an additive bias larger than ^ 10^"^ would degrade the 
cosmological constraints of future weak lensing probes. This re- 
sult suggests that whenever the PSF is not well corrected and the 
residual is constant aperture mass statistics should be preferred. 

We show in the right panel of Figure|4]a more realistic case of 
modified top-hat variance, which assumes the PSF is modeled in- 
dependently in sub-regions of the image. In this case, the additive 
bias only correlates at the scales of those sub-regions. In this exam- 
ple we assume the correlation is constant below ~ 5', which is a 
typical scale of the CCDs of ground-based cameras, and decreases 
steeply (exponentially) to zero at larger scales. In this toy-model 
the PSF residuals are strongly correlated at small scales enhancing 
strongly the signal, whereas on large scales > 5' the PSF resid- 
uals are decorrelated and only the multiplicative constant affects 
the signal. We expect something similar to happen on real data as 
the PSF varies across the field and it is generally corrected by in- 
terpolating the PSF shape over a characteristic area. If this is the 
case, the additive constant would affect both the top-hat and the 
aperture mass variances and the effect could be significant at small 



scales. However, the STEP2 simulations have a PSF which is con- 
stant across the field of view and this does not allow one to explore 
in more detail the effect of such a residual on real data. 



4 PARAMETER SPACE DEFINITION 

In order to estimate the impact of the redshift-dependent shear mea- 
surement bias on cosmological constraints, we perform a likelihood 
analysis in a nine-dimensional grid using two-point shear statistics 
with and without tomography. 

We define the log-likelihood of a model as 

C{ch) = i[d.,W-«.,(tf;0)rC(^,t?';<^)-i 



(11) 



Both the tomographic vector of data, dij(j}), and the two-point 
functions theoretical predictions, Vij {•&; <j)), are computed at twenty 
angular scales, between 6 "and 1 °, and correlate three redshift 
bins i,j = [1, 3]. The covariance matrix of the estimator used to 
measure the two-point shear statistics, C(i9, (j>), has dimension 
120 X 120. 

The redshift bins are the ones used in Fig. |4l = 
[0,0.77], Z2 = [0.77,1.15], 23 = [1.15, 4.0], with mean red- 
shifts of zi — 0.47, Z2 — 0.94, 23 = 1.69, respectively. They 
were chosen such as to optimise the information. Since the sig- 
nal increases with redshift, it is more efficient to divide the red- 
shift range such that the lower-redshift bins contain a larger frac- 
tion of the galaxies. We followed the strategy of Hu (1999) and 
defined the bins such that 21 contains half of the galaxies and 22 
and 23 contain one quarter each. The correlation coefficient be- 
tween the aperture mass variances of the various bins is defined 
asn,(^) = M2p,,,(t9)/(Mfp,„(i?)Mfp,^^(i?))i/2. We measure 
a large correlation ri2 ~ r23 ~ 0.9 between consecutive bins, 
which shows that further subdivisions do not increase the amount 
of information. 

The data vector, dijij)), is computed for a fiducial model 
similar to the WMAP5 average values (Q.m ~ 0.27, Qx — 
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0.73, h = 0.705, 0-8 = 0.812, ^6 = 0) for a ACDM cosmol- 
ogy ( [DunMey et al. 2008) and no bias. 

The value of Vij[(j>) depends on cosmological parameters 
and also on the bias parameters a^, fem, Oc, b^. The chosen range 
of the bias parameters covers the results of the STEP2 analysis: 
am = [-0.03,0], &™ = [-0.06,0.06], = [-0.001,0] and 
be — [—0.006, 0.002]. The cosmological parameters varied in the 
grid are: the matter density Q^n ~ [0.1, 1.0], the Hubble constant 
h = [0.6, 0.8], the normalisation of the matter power spectrum 
(78 = [0.5, 1.1], and two parameters, wo and wi, for the equation 
of state of dark energy, u)de(z) = Pve{z) / Pde{z) ■ We model 
wdb(2) in two different ways. The first one is a constant equa- 
tion of state, woe{z) = wq = const; for this model we vary 
Wo = [—2.0,0]. The second model has been suggested by Ben- 
abed & van Waerbeke (2004): 

C ) - f ^0 +wi log(l + z) , z ^ 1 

\ Wo + TOi[Iog(2) — arctan(l) + arctan(2;)] , z > 1 

(12) 

Benabed & van Waerbeke (2004) have shown this model is able 
to mimic quite accurately the Universe dynamics generated by 
SUGRA potentials. For this case we consider the ranges wo = 
[-2.0,0], wi = [0,0.4]. 

For the source redshift distribution we use the parameteriza- 
tion of Eq. ([9j. We do not vary its parameters in the likelihood anal- 
ysis. This choice is supported by the fact that the errors affecting the 
estimation of the parameters a, (3 and zo are about one percent, in- 
cluding the sample variance affecting the CFHTLS-Deep catalogue 
( Ivan Waerbeke et al. 2006J . 



5 COVARIANCE MATRIX 

In this section we compute the covariance matrix C(i?, (j>) for 
Eq. din . Most studies of weak lensing tomography have been made 
in Fourier space where the covariance matrix is easy to compute. 
Analyses made in real space usually evaluate the covariance ma- 
trix from Monte Carlo simulations (eg., Simon, King & Schnei- 
der 2004), or more recently analytically from the covariances of 
power spectra jJoachimi, Schneider & Eifler 2008^ . Here we fol- 
low a different approach deriving analytically in real space C++ = 
C(^+,ij('i}i),S,+.ki{'&2)), the covariance matrix of 5+,ij(i9), ex- 
tending to tomography the formulae derived in Schneider et al. 
(2002). The details of the derivation are given in appendix A. The 
covariances of the estimators of {A/fp(i9)) and {7^(1?)) are af- 
terwards derived from the covariance of C+(i?) (Schneider et al. 
2002). 

The pure cosmic variance term is analytically derived in Ap- 
pendix A assuming a Gaussian shear field [Eqs. l |A9t and dAlOU , 
for the a fiducial model defined in this case defined in the previous 
section. We then recalibrate this matrix for non-Gaussian contribu- 
tions following Semboloni et al. (2007), hereafter S07. This con- 
sists in multiplying the Gaussian cosmic variance by a calibration 
factor ij; i?2, kl). 

We found that the calibration factor found in S07 needed to 
be modified for the case of tomography by applying it to Eqs. l |A9| l 
and l lAlOb . and comparing the result with the non-Gaussian cosmic 
variance for tomography, measured directly on the same set of ray- 
tracing simulations used in S07. The best agreement is obtained 
when the redshift used in the calibration factor of S07, 

^('^^'^r,^2,ki)^^^^^^, (13) 



is redefined as an effective redshift z — z{ij; kl) given by 

z{ij; kl) = z= \ \J{zi + Zj) {zk + zi) . (14) 

The effective redshift is thus the geometric mean of the aver- 
age redshifts of the bins involved in each of the two correla- 
tions ^+,ij and 5+,fci. The functional form of a{z) and [3{z) 
are defined in S07, with the same best-fit parameter values 
{a\,a2, as; bi, 62, 63, 64; ti, i2, is) indicated there. Finally, the re- 
calibration is only applied at angular scales below the J?max defined 
in S07. 

We find that this fitting formula agrees with the measured co- 
variance to an accuracy of 20%; the worst result being for the lower 
redshift bin. This recalibration might overestimate the values of 
the covariance matrix l ITakada & Jain 2008} (as the simulations we 
used have a high value as, — 1) but this is not so important as we 
include in our likelihood analysis large scales measurements. 

After recalibrating the cosmic variance term, we add the other 
contributions to C++: the coupling between shape-noise and the 
shear signal given by Eq. l lA12t . and the pure shape-noise term 
given by Eq. l |A13| l. These terms depend on the density of galaxies, 
the total area of the survey and the intrinsic ellipticity dispersion. 
We assume a density of galaxies integrated along the line of sight 
of n = 12/arcmin'^ which is similar to the CFHTLS-Wide sur- 
vey with a median redshift of Zm = 0.8. To find the density of 
galaxies for each redshift bin, we normalise the redshift distribu- 
tion of the sources, p[z), to n and integrate it between each ^int 
and Zsup. We use an intrinsic ellipticity dispersion CTj = 0.44. The 
cosmological constraints are derived for two different survey areas: 
A — 170 deg'^ and A = 2000 deg^, which are a good represen- 
tation of current and future weak lensing ground-based multicolor 
surveys. The first case has the same sky coverage and depth as the 
CFHTLS-Wide survey. The second has similar sky coverage to the 
upcoming KIDS survey, while being about half a magnitude deeper 
than KIDS. However, as the shape measurement bias is a function 
of signal-to-noise, we expect it to be larger, for the same magnitude, 
in a shallower survey than in a deeper one. 

Finally, using Eq. (42) of Schneider et al. (2002) with K+ — 
1, we obtain the covariance matrix of the estimator of {A'/fp(i9)) 
from C++, and using a similar expression we derive the covariance 
matrix of the estimator of (7^(1?)). 

During the likelihood analysis the covariance matrix is kept 
constant, neglecting its dependence on the cosmological parame- 
ters. The accuracy of this approximation is studied in Eifler, Schnei- 
der & Hartlap (2008). We also neglect its dependence on bias pa- 
rameters. Thus in the likelihood defined by Eq. i ll It . the depen- 
dence on the bias is only contained in Vij . This type of systematic 
error which affects the signal but does not introduce extra noise in 
the covariance matrix has been called Type I by Kitching, Taylor & 
Heavens 2008a and it allows for self-calibration. 



6 RESULTS 

In this section we obtain cosmological constraints by performing 
the likelihood analyses previously described, both with and without 
tomography. As stated earlier on, our approach is to include the 
realistic bias measured in STEP2, with as little extra modeling as 
possible. For this reason we choose not to model the correlation of 
the additive term, which does not come directly from the STEP2 
analysis. We will focus thus on the multiplicative bias and show 
constraints using the aperture-mass dispersion only. 
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Figure 5. Cosmological constraints in tiie [Sim, erg] plane for a survey of area 170 dcg-^ and three tomograpliic bins. The contours represent the 68% and 95% 
confidence regions and are marginalized over the hidden cosmological pai'ameters. Left panel: the effect of a fixed bias of am = —0.02, bm, = (dashed, 
blue) is compared with the bias-free contour (sohd, red). Right panel: the biased contour is now marginalized over the bias parameters (dashed, blue) and 
compared with the bias-free contour (solid, red). 
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Figure 6. Confidence regions (68% and 95%) on the [Qm,wo] and 
[Qm, o-g] planes for the 2000 deg^ survey. The contours, marginalised over 
the hidden parameters, are shown for the case of no bias (solid, red) and 
marginalising over the shape measurement bias parameters (dashed, blue). 
Left panels show results of the non-tomographic analysis and right panels 
refer to the tomographic case. 



We start by considering the 170 deg survey, which has rela- 
tively large noise, to asses if a redshift-dependent shape measure- 
ment bias already affects the results of a tomographic analysis of a 
survey such as the CFHTLS-Wide. The constraints in the [fl 
plane, using auto- and cross-correlation in the three redshift bins 
previously described, are shown in Figure|5] In the left panel, con- 



tours are shown for the case of no bias and a particular choice of 
bias {am = —0.02 , bm = 0). The biased contour is clearly shifted 
to the upper-right comer to compensate the negative multiplicative 
bias. This shows that our choice of redshift binning is quite effec- 
tive, having a noise level low enough for systematics to become 
important. The shift in the constraints can be quantified by fitting 
the direction of the degeneracy seen in the plot. In the case of no 
bias the fit is given by 

(78(^^/0.275)°-'** = 0.79 ± 0.02 , (15) 

whereas in the case of this particular bias we find 



a8(J7m/0.275)° 



0.81 ±0.02. 



(16) 



This result shows that ignoring the presence of the shape measure- 
ment bias would lead to biased contours. 

When the amplitude of the bias is not known, the safer 
approach is to marginalize over a reasonable interval of bias. 
Marginalizing over the range suggested by our STEP2 analysis, 
we find the contour shown in Figure |5] (right panel). The preci- 
sion of the tomographic measure of [Tim, cs] is reduced by a factor 
of 2. This result includes the degeneracies between cosmological 
and bias parameters. Marginalizing over the cosmological param- 
eters, the bias parameters could be estimated from the data, i.e., 
the signal could be self-calibrated. For the self-calibration to be ef- 
fective, priors would be needed. They can be obtained with extra 
information from independent measurements affected by the same 
systematics. Other measures of cosmic shear such as higher-order 
correlation functions ('Bemardeau, van Waerbek e & Mellier 1997| l 
or the shear-ratio test (Jain & Taylor 2003 1 could be helpful for this 
purpose. 

Finally, we also produced contours in the [iyo,r2m] plane. 
In this case the systematics do not have a relevant effect, since 
the bias-free contours are too broad and the estimate is noise- 
dominated. 

We turn now to the case of the survey of A = 2000 deg^, in 
order to see how the shape measurement bias can affect the weak 
lensing analysis of data from future ground-based observations. We 
quantify the constraints using a figure-of-merit, FoM, defined as the 
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Figure 7. Marginalised 1 -dimensional probability distribution of wq for 
two cases: no bias (solid line, red) and a shape measurement bias character- 
ized by am = -0.02, 6™ = (dashed hne, blue), for the A=2000 deg^ 
survey. The peak of the probability distribution changes significantly with 
the shape measurement bias parameters. 

inverse of the area enclosed by the 68% confidence level. Quoted 
results imply a marginalisation of the bias over its full range. 

In Figure |6l we show the cosmological constraints in the 
[nnDCs] and [f2ni,™o] planes, both for tomography and non- 
tomography, using the variance of the aperture mass. The first in- 
teresting result is that tomography produces an enormous gain in 
all constraints, but especially in dark energy. Looking at the bias- 
free contours, the increase in the FoM between using and not using 
tomography is a factor of 10 for the [S7m,TOo] contour and a fac- 
tor of 6 for [Qm, os]. This is larger than the improvement between 
the two surveys due to the increase in area, which scales with the 
square root of the survey area and is thus around 3.5. 

Regarding the impact of the bias, we find the multiplica- 
tive redshift-dependent bias essentially does not affect the non- 
tomographic constraints, even for a KIDS-like survey. Indeed, the 
FoM of the [f2m,wo] contour only decreases by a factor of 1.1 
in the presence of bias. For the [fimjCrg] contour the decrease is 
larger, around 1.6. The tomographic constraints are however are 
however more strongly affected by the bias. We find using tomog- 
raphy the [Qm, erg] marginalised FoM is ~ 3.5 smaller in compar- 
ison to the tomographic FoM without bias. It is still better than the 
bias-affected non-tomographic one by a factor of around 2, but it is 
only slightly better than to FoM for the A = 170 deg^ survey in 
the case of no-bias. Notice however that even though the two FoM 
are similar, the corresponding one-dimensional constraints for Qm 
and as are different. Indeed the two contours (the solid contour in 
the right panel of Fig.|5]and the dashed one in the upper-right panel 
of Fig.|6f have different shapes. This indicates the effect of the bias 
is not equivalent to an increase of the data error bars. Similarly, 
marginalising over the bias causes the FoM of the tomographic 
[f2m, Wo] constraint to increase by a factor of 2 in comparison to 
the unbiased FoM. 

The finding that tomography is more affected by the bias re- 
flects the fact that cosmological and bias parameters are more cor- 
related for this case. We looked at the bias-cosmology correla- 
tions and verified that while for the non-tomographic case only Q,n-i 



shows some correlation with am and bm, in the case of tomogra- 
phy Sim, erg and wo are all correlated with am - In particular, erg and 
a,„ are anti-correlated, whereas wo and a,„ correlate. This means a 
lower (more negative) wo is needed to compensate the multiplica- 
tive bias, which essentially underestimates the signal. This result 
can be explained if one considers how the properties of large-scale 
structures depends on dark energy. For larger wq, dark energy starts 
to dominate at times later than those predicted by ACDM, and the 
structures grow faster. If one normalises the amplitude of the den- 
sity fluctuations to the present time then models with uiq > — 1 
would have less structure in the past and produce a signal ( A/fp {■&)) 
which is smaller than the signal produced by the ACDM model. 
The dark energy affects two-point cosmic shear statistics also via 
the geometric factor and this effect has the opposite result increas- 
ing the amplitude of {Mfp(i9)) with wo ( [Simpson & Bridle 2065] l. 
This last effect is dominant for the correlations involving the lower 
redshift bins, whereas the previous effect is dominant for high red- 
shift bins. For our particular choice of redshift bins, the first effect 
is the dominant one, as is shown in Figure|7] by the shift to the left 
of the probability distribution, in the presence of bias. This plot also 
illustrates that ignoring the presence of a shape measurement bias 
on the likelihood analysis may lead to an incorrect cosmological 
interpretation. 

We also performed a likelihood analysis for the A — 
2000 deg'^ survey, this time using the equation of state of dark 
energy given by Eq. l ll2t . including an extra parameter, wi. The 
resulting marginalised contours for both the [f2m,(Tg] and the 
[flm, Wo] plane are very similar to the results of the previous anal- 
ysis. The left panel of Figure [8] shows the contours in the [too, wi] 
plane. They are weakly affected by the bias, broadening mildly 
towards lower values of wo consistently with the previous anal- 
ysis. {A/fp(i9)) cannot independently measure wi, as shown by 
the rather flat probability distribution in the right panel of Figure[8] 
(solid line). However for a bias of {am ~ —0.02 , bm ~ 0), i.e., 
assuming this bias in the theoretical models and keeping the fidu- 
cial model bias-free, the probability distribution is much narrower. 
The reason for this is that the impact of this bias on the models is 
much larger than the effect of wi. Therefore, most values of wi 
cannot compensate for the bias impact. Conversely, when doing a 
likelihood analysis of real data biased by an unknown amount, as- 
suming no bias in the theoretical models, the same effect may occur 
leading to an inaccurate but apparently precise estimate. 



7 CONCLUSIONS 

We have investigated the impact of realistic redshift-dependent 
measurement bias on the estimation of cosmological parameters 
using two-point shear statistics in a tomographic approach with 
three redshift bins. We focused on two survey types of the same 
depth, with median redshift ~ 0.8, covering an area of 170 deg^ 
and 2000 deg^, representing current and near-future ground-based 
weak lensing surveys, respectively. 

We performed a likelihood analysis in a grid of cosmologi- 
cal and bias parameters. The covariance matrix of the tomographic 
two-point functions was derived analytically in the real space, ex- 
tending the formulae of Schneider et al. (2002), and was after- 
wards calibrated for non-Gaussianity, extending the fitting formula 
of Semboloni et al. (2007). 

Realistic redshift-dependent multiplicative bias were obtained 
reanalyzing the STEP2 set of simulations with a KSB pipeline. The 
results of our analysis are thus strictly valid for KSB -based mea- 
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Figure 8. Cosmological constraints on the parameters of the dark energy equation of state, for the survey of area A = 2000 deg^. Left panel: 68% and 95% 
confidence regions on the [wg, wi] plane for the non-bias case (solid, red) and marginalised over the bias parameter space (dashed, blue). Right panel: the 
marginal 1-dimensional probability distribution for wi obtained assuming no bias (solid line, red) and assuming a fixed bias of {am = —0.02 , = 0) 
(dashed line, blue). 



siirements only. In this framework, our results are also a worst-case 
scenario, since the higher end of the bias values found correspond 
to the highly elliptical PSFs D and E of STEP2, which will be 
avoided by design in future surveys. Nevertheless, such elliptici- 
ties are comparable with those measured in the outer regions of 
CFHTLS-Wide fields. 

The top-hat variance is more sensitive to the presence of a 
shape measurement bias than the variance of the aperture mass. In 
particular, in the case of a large scale-independent additive bias, 
cosmological interpretation of the signal may even not be possible, 
suggesting that aperture mass statistics should be preferred to top- 
hat statistics. The measurement of {Affp(i?)) is however also quite 
sensitive to the shape measurement bias; its amplitude is reduced 
by 20% at redshift ~ 1, and its scaling with redshift becomes shal- 
lower. 

We stress the results shown in this paper are based only on 
the simulations. Using these simulations the aperture mass variance 
seems to be more accurate than the top-hat variance, this because 
the top- variance is particularly sensitive to the additive bias as long 
as this bias is constant across the field. On real data this is not often 
true and we expect both two-point shear statistics to be biased. Un- 
fortunately, we could not quantify this effect using this set of sim- 
ulations as the PSF is constant across the field. Moreover, in real 
data, the profile of the PSF is not constant in time, space and color, 
and it might have a more complex shape. These factors may intro- 
duce other sources of bias which are not investigated here. Finally, 
in real data, the lack of accuracy on the estimation of photometric 
redshifts could further decrease the scaling of the two-point shear 
statistics with redshift. Indeed, the presence of outliers in photo- 
metric redshift estimation is known to decrease the measured slope 
of the evolution of the cosmic shear signal dErben et al. 2008l l. 

We found the redshift-dependent bias to have, in general, a 
large impact on the tomographic measurements. For example, con- 
straints obtained from the largest survey, marginalised over the re- 
alistic bias interval, are comparable to the ones from the ten times 
smaller survey assuming the latter does not suffer from measure- 
ment bias. 

Ignoring the presence of a shape measurement bias can al- 
ready bias the estimation of cosmological parameters in surveys 



such as the CFHTLS-Wide. From our analysis, the bias on erg is of 
a few percent. For the survey of 2000 deg'^ the effect of ignoring 
the bias is more important, around 5%, corresponding roughly to 
1 a, for Wo, and larger for wi. 

The safer approach when dealing with unknown biases is to 
marginalise over a realistic range. While this reduces the constrain- 
ing power of the tomographic analysis over a non-tomographic one, 
there is the possibility for self-calibrating by combining informa- 
tion from other cosmic shear measures. However, the best approach 
would be to correct for these biases. For such a reason, collabora- 
tions that aim to improve the measurement of the galaxy shapes 
such as STEP and GREAT08 dBridle et al. 2008l l are important for 
the development and improvement of shear measurement methods. 
As an example, the simulations produced by the STEP collabora- 
tion have been used to develop and test new PSF-correction meth- 
ods, such as shapelets-based decompositions ( |Kuijken 2006 |l and 
the LENSFIT method ( [Kitching et al. 2008b| l, which are among the 
most promising methods for future weak lensing analysing. 

To conclude, this analysis, which is limitated to the KSB PSF 
correction method and uses only the STEP2 set of simulations 
shows the importance to investigate further the dependence of the 
shear measurement bias on the simulated galaxy properties in order 
to infer more accurate constraints on cosmological parameters. Fur- 
thermore, simulations cannot be used to investigate all the sources 
of systematic errors which may be present on real data; for this 
reason one needs to complement this study using consistency tests 
with real data, for example measuring the cosmic shear signal of 
various galaxy populations with different color, size, magnitude 
and ellipticity. 
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APPENDIX A: COVARIANCE OF THE TOMOGRAPHIC ESTIMATOR OF ^4 



We extend the calculations of Schneider et al. (2002) (S02 hereafter) in order to explicitly include the redshift binning of the correlation 
function. Note latin subscripts {i, j) refer to individual galaxies, while greek subscripts (a, j3) refer to bins. 

The shear 7 at angular position Oi of a source galaxy at redshift Zi is estimated from the observed ellipticity of the galaxy image. After 
deconvolving the PSF, the corrected ellipticity is related to the intrinsic ellipticity and the shear by 



ei ^ el + jiOi, Zi) 



(Al) 



The correlation function Z{a, [3)) between redshift bins with source redshift distributions Pa{z) and P0{z) is estimated in 

angular- separation bins The angular bins are defined by the function A^{\9i — 9j\) ~ 1 for 1? — Ai?/2 < \9i — 9j \ ^ ■& + A-3/2 and 
zero otherwise, while the redshift bins are defined by Az(zi, Zj) = 1 for pa{zi) x {zj ) 7^ and zero otherwise. In other words, only pairs 
of galaxies with one of the galaxies belonging to the bin a and the other to the bin /3 contribute to the correlation at Z{a, 13). The estimator 
used for the correlation function is. 



' Y.-^WiWj {ei+ej+ + eixejx) A^QOi - 9j\)Azizi,Zj) 

^^^"'^^ = Ju^) 



(A2) 



where the subscripts (+, x) refer to the two components of the ellipticity, Np{-d, Z) is the effective number of pairs contributing to the 
correlation at Z) and Wi, Wj allow for weighting the galaxies. In the absence of intrinsic alignments and shear-shape correlations, 
Eq. iA2l is an unbiased estimator of The noise of the estimator is a function of the ellipticity shape noise, (ei+ej^ + e^x ^jx ) — f^^ij ^ 
which we assume is independent of redshift. 
The covariance of ^+(1?, Z) is defined as 

C0Y{i+ ,d^,Zv,i+,^2,Z2)^{i+ (^1 , ^1 ) i+ (^2 , ^2 )) - C+ , 2l ) C+ (^2 ,Z2). 

Inserting Eq. l lA2t it reads, 

Cov(|+, j9i,2i;|+,i92, Z2) = 'Y^WiWjWkWiA^^{ij)Ao^{kl)Azi{ij)Az2{kl) {{(-ntji + ei2ej2) (efciQi + efc2Q2)> X 

ij k I 

1 



(A3) 



iVp(l9l,2l)7Vp(^2,22) 



-5+(^l,2l)C+(^2,22 



(A4) 



This expression depends on four-point correlations involving the four positions {9i, Zi), {8j, Zj), {9k, Zk) and {61, zi), where the first two 
positions define (i9i,.Ei) and the latter two define (192,^2). In order to evaluate it, we proceed as in S02: inserting Eq. l lAlb : assuming a 
Gaussian shear field and thus factoring four-point functions as a sum over products of two-point functions (assuming no intrinsic alignments), 
and finally writing the two-point correlations of the shear components as function of ^+ and f _ . We obtain the results (A5)-(A7), which differ 
from Eq. (23) of S02 since we can no longer interchange i with j or k with which may have different redshifts. The cosmic variance term 
reads 



V++(e+,l;^+,2) 



Np(^i,Zi)Np{^2,Z2) 2 



WjWkWiA^^ (zj) A^2 (^O^^i {'^3)^2,2 (^0 



ijkl 



i+{il)^+{jk) + 



+ i+(ik)£,+ {jl) + cos[A{ipa - ip,k)]i-{il)i-(jk) + cos[A(ip,k - 1/9^;)] C- (ifc)C- O'O 



(A5) 



where tpij is the polar angle of the vector 9i — 9j. The mixed term, with the coupling between shape-noise and shear, is now 



M++(C+,l;e+,2) = 



^ wiwjWkAi,^ (y ) (ik)Az^ {ij)Az2 {ik)^+{jk) + 



Np{i)i,Zi)Np{^2,Z2) 2 

+ ^ w.WjWkAi,^ {ij)^«2 {jk)Azi {ij)Az2 ijk)^+{ik) + ^ wfwjWiA^f^ {ij)^«2 {il)^Zi {ij)Az2 (iOC+O'O + 



ijk 



ijl 



+ ^ w^WJWlAi,-^ {ij)Ai)2 {jl)Azi {ij)Az2 (iOC+(»0 

ijl 



(A6) 



The noise term is similar to the non-tomographic case, only contributing to the diagonals of the covariance matrix, i.e. for i9i = ^2 and 
Zi — Z2. It reads. 



S(C+,l;e+,2) = 



Np{i}i,Zi)Np{i}2,Z2 



■at ^ WiWjA^)^ (ij) Ai>2 (i3)Az^ iW^Zi iW 



(A7) 



Table |A] summarizes the dependence of the cosmic variance on the tomographic auto- and cross-correlation functions, considering two 
redshift bins for illustration purposes. In the case the cross-correlation signal is not considered in the data vector, the matrix reduces to the 
sub-matrix highlighted, which nonetheless depends on the tomographic cross-correlation. Note this restricted case corresponds to imposing 
Pa{z) = Pi3{z) being thus invariant under i — j and k — I transformations, and therefore Eqs. JASt -l lATt reduce to Eq. (23) of S02. 
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€ii C22 C12 



?11 


+s 


?12 5l2 


£11 £12 


?22 


5l2 C12 


?22 ?22 + S 


£22 £12 


?12 


£11 £12 


?22 £.12 


£11 £22 + £12 £12 + s 



Table Al. Summary of the tomographic cosmic variance matrix, with dimension 3x3 con'esponding to 2 redshift bins {a = [1, 2]), showing the dependence 
on the various correlation functions £-^-{"3, Z{a, (3)) (denoted by £cp)- The terms affected by shot noise are also indicated (with S). Dashed lines isolate the 
restricted case of neglecting the information in the cross-correlation of redshift bins. 



The calculations in Eqs. JA5t -l lA7] l are time-consuming and survey-dependent since they involve sums over the galaxies positions. Still 
following S02, we calculate a more convenient quantity: the ensemble average of the covariance matrix for a survey of area A and galaxy 
density n, and write the result in the form 

E (^Cov(|+,t?i,2i(ai,/3i);C+,i?2,^2(a2,/32))) = D So{^i ~ ■&2) SuiZi - Z2) + q++ + r+o + r+i. (A8) 

For the cosmic variance we obtain 
1 



r+o = 
r+i = 
+ 



■kA 



d(t)(t) 



dvPlC+,QiQ2(l'0al) / (, + , 13^02 b\) + dvSlC+.aifedV'al) / d<^2 C+ ,/3i ( I '/'b I ) 



(4^)^ Jo 

1 



d(t)(j) d(y3l 5_,QiQ2(l'0al) / dv32C-,/3lfe(l^bl) [C0s4v5a C0S4((96 +Sin4v5a Sin4((96] 



(A9) 
(AlO) 



(4^)^ Jo 



d(l)(j) d(^l ^-.Qi^gad'/'al) / d((92C-,/3ia2(|^6l) [C0s4v5a COs4v3b +Sin4v3a Sin4((96] 



where ^+,ap is a short-hand notation for ^+{^, Z{a, (3)). The result was obtained as in S02: the integrations arise from the averaging 
operator, and a change of variable is made from the position angles 6, included in Eq. iA5i , to separation angles (j). We note that due to the 
loss of the i — j, k — I invariances, different changes of variables are needed for the two terms of r+o or r+i ; (f) is defined (f) — 9s ~ 6i 
in one case and (j) = 63 — 62 in the other. The angular separation information is contained in the vectors l/j^ and their polar angles ipa 
and ipt'. 



<j)cos ip ~ "di cos ifii 
(f) sin 99 — t?i sin 991 



4>COS ((5 + 1?2 cos ((32 

(f) sin ^ + 1^2 sin ip2 



(All) 



For the mixed term we find 



1 

2-kA 



d</5 



:,/3l/32 



^+,c«ic2 2l?l1?2 COS if] + 



V"'/3l"/32 



e+,/3l/32 J ^1+ ^1-2^1^2 COS cp + 



2 

""e,/3ia2 c 



The expectation value of the shot noise term is 

D = 



2 + i?i - 2did2 cos <^ + 



A47ri9 Ai9 



""e,/3i/32 



1?? + 1?^ - 2l?i1?2 COS<y3 



+ 



(A12) 



(A13) 



2 2 2 



where Ai9 is the angular bin size. We assume o-^ ^^^j.^ = '^'t,f3ifi2 

Finally, we insert Eq. iA8t in Eq. (42) of S02 to obtain the covariance matrices for the aperture-mass and top-hat dispersion. In the 
absence of B-modes they depend only on the covariance of hence we do not calculate Cov(^_ \^~) or Cov(^+; f 



